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Abstract. In order to investigate corrections to the common KdV approxi- 
mation to long waves, we derive modulation equations for the evolution of long 
wavelength initial data for a Boussinesq equation. The equations governing 
the corrections to the KdV approximation are explicitly solvable and we prove 
estimates showing that they do indeed give a significantly better approxima- 
tion than the KdV equation alone. We also present the results of numerical 
experiments which show that the error estimates we derive are essentially op- 
timal. 



1. Introduction 

Modulation, or amplitude, equations are approximate, often explicitly solvable, 
model equations derived — usually through asymptotic analysis and the method of 
multiple time scales — to model more complicated physical situations. Although 
these equations have been been used for over a century, only lately has there been 
an attempt to rigorously relate solutions of the modulation equations to the original 
physical problem. In particular, through the work of Craig m, Kano and Nishida 
13 1, Kalyakin Schneider Ben Youssef and Colin |^ and Schneider and 



Wayne , |^ | , the validity of Korteweg-de Vries (KdV) equations as a leading 
order approximation to the evolution of long wavelength water waves and to a 
number of other dispersive partial differential equations has been established. 

While the KdV approximation is extremely useful due to its simplicity and the 
fact that the KdV equation can be explicitly solved by the inverse scattering trans- 
form, both experimentally and numerically one observes departures from the predic- 
tions of the KdV equation. Our goal in this paper is to derive modulation equations 
which govern corrections to the KdV model. In the present paper we will not work 
with the full water wave problem but rather study modulation equations for long 
wavelength solutions of the Boussinesq equation: 

dtt ~ dxx — {d'^)xx + dttxx, 

a; e K, i > 0, e{x,t) e R. 

Our motivation for studying this equation is twofold. First, the Boussinesq 
equation was originally derived as a model equation for water waves, and as our 
ultimate goal is to derive corrections to the KdV approximation to water waves, 
we regard the study of (|l|) as a useful first step in understanding the much more 
complicated water wave situation. We note that Schneider's analysis of the KdV 
approximation for (|l|) in pl| ] served as a template for the analysis of the water wave 
problem in p2[. 
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Our second justification for deriving second order modulation equations for ([^) 
is that these modulation equations serve as a sort of normal form for more compli- 
cated PDEs, and as such we expect that the modulation equations which describe 
corrections to the KdV approximation to ([|) will also govern corrections to the 
KdV approximation in more complicated situations as well. Thus the results on 
existence, uniqueness and other properties of the modulation equations we derive 
in this paper should also be of use in more complicated situations that we plan to 
treat in the future. 

We now describe in more detail our results. It is convenient to rewrite as a 
system of two first order equations. As in |2l| we introduce new variables, 

u{x,t) ^l{e{x,t) - X-^0t{x,t)) 
(2) I 

v{x,t) ^-{e{x,t)+X-^9t{x,t)) 

where A is a skew-symmetric multiplication operator in Fourier space defined by 
Xu = (ifc/Vl + Note that for X~^9t to be well-defined, we must have 9t{0, t) = 

0. That is, the average value of 9t should be zero. We note that for 9{x,t) a solution 
of (0), we have that the average value of 9t{x, t) is a constant of the motion. Thus, 
if the initial condition for 9t has zero average, 0t(O,t) = will remain zero for all 
time. Furthermore, as discussed in assuming that the initial condition 0t(x,O) 
has zero average is not unnatural considering the origin of (0). Thus, we will make 
that assumption so that the change of variables (||) is well-defined. 
Taking time derivatives of u and v we find: 

u \ f —A \ f u \ If —X{u 



^'{v J~{ X J{v J^2{ X{u + vr 

Not only is (|3[) convenient from a mathematical point of view but, as we shall 
see, u and v have the physical interpretation of being the left and right moving 
parts of the solution. 

We turn now to the assumptions on the initial conditions of (||). The KdV 
equation is an approximation of small amplitude and long wavelength motions, 
and thus we will assume that the initial conditions of (|l|) are of this form. More 
precisely, fix a constant C/ > and assume: 

Hypothesis 1. There exist Uo, Vo with 

max{|lf7o||//=(4)n//=+f, l|Vb||ff=(4)nff=+9} < Ci 
such that the initial conditions of are of the form: 

u{x,0) = e^t7o(ex) , v{x,0) = €^Vo{ex) , 

where e is small. 

Here, iJ^(m) = {/|(1 + a;^)"/^/ e H'}. The norm on this weighted Sobolev 
space is given by 1 1 / 1 1 jj. („) = 1 1 ( 1 + ) / 1 1 ■ We use these spaces because we are 
interested in solutions which are in some (weak) sense "localized" . In particular, 
any small perturbation of the known soliton solutions to the KdV equation will 
satisfy this localization property. 

Remark 1. We can, of course, recover the initial conditions for ^) from Hypoth- 
esis ^ via and we see that the initial conditions expressed in the 9 variables are 
also of small amplitude, long-wavelength form. 
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According to the KdV approximation results of long wavelength solutions 
of split up into two pieces, one a right moving wave train and one a left moving 
wave train. Each of these wave trains evolves according to a KdV equation and 
there is no interaction between the left and right moving pieces. One might expect 
two types of corrections to such an approximation: 

• corrections due the fact that the left and right moving wave trains will interact 
at higher order. 

• corrections due to the fact that even in the case of a purely right (or left) mov- 
ing wave train, solutions to the Boussinesq equation are not exactly described 
by solutions to the KdV equation. 

Both of these types of corrections are apparent in our results and in fact the correc- 
tions to the KdV approximation are a sum of solutions of two types of modulation 
equations; an inhomogeneous transport equation and a linearized KdV equation 
which can be seen (roughly speaking) as reflecting these two sources of corrections. 

To incorporate these two types of corrections, we add to the KdV wavetrains, 
which we denote by U and V (since they represent the leading terms in u and v 
respectively), additional functions A and B and F and G. These functions then 
satsisfy the modulation equations: 

(4) ,2 2 

dTV= -d\^V+-dxX 

drA + dxA = - ]-dx^ V^iX + T, eV) - dxU{X - t, eV)F(X + r, eV) 

(5) { 

drB - dxB ^ -dx_ U^{X - r, eV) + dxU{X ~ r, eV)F(X + r, eV) 

and 

OtF = - dx_ (UF) - iai F + ji 

(6) ,2 

5tG= dx+{VG) + -d'x+G + 

where T — e't, t = er, X — ex and X± = X ±t. 

The first of these pairs of equations is simply the KdV approximation. The 
second and third pairs give rise to the corrections to the KdV approximation. We 
note that the terms and which appear in (^) are inhomogeneous terms which 
are made up of a combination of sums and products of the solutions to (^, (||) and 
their derivatives (see equations (p^-(|l6|) below). 

There is some freedom in how we choose the initial data for the modulation 
equations. For simplicity we assume that U{X,0) = Uo{X), V{X,0) — Vo{X) 
and choose zero initial data for (|) and (|), i.e. A{X,0) = B{X,Q) = F{X,{)) = 
G(X,0) = 0. 

That the KdV equation has solutions for all times with initial data of the type 
described is well known. In particular one has (see [p2[): 

Theorem 1. Let s > A. Then for all Cq,Tq > there exists Ci > such that if 
U , V satisfy (0j with initial conditions Uq, Vq and 

max{||C/olU=(4)nff=+9, || Vo|U=(4)nff=+9} < 
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then 

(7) sup {||C/(-,T)||H=(4)nH=+«,im-,r)||H.(4)nH^+s} < Ci 

Te[0,To] 

On the other hand it is less clear that solutions of (|^) and (||) will remain bounded 
over the very long time scales necessary for the KdV approximation. Thus, the first 
significant technical result of this paper is: 

Proposition 1. Fix To > 0. Suppose, C/o, Vb G F'^(4) and U, V, A, B, F and G 

satisfy then there exists a constant C2 such that the solutions of and 

^ satisfy the following estimates: 

sup {||A(.,r)||ff.-3,||B(.,r)||H.-3} < C2 

rG[0,To£-2] 

sup {||F(.,r)||^,||G(.,r)||^} < C2 

Te[0,To] 

where H = H'^-^ D H''-^{2). 

With this preliminary result in hand, we can now state our principal result. 

Theorem 2. Fix To, C/ > 0, a > 13. Suppose U, V, A, B, F and G satisfy 
equations ^)-^- Then there exists eo > and Gp > such that if the initial 
conditions for satisfy Hypothesis^ then for e G (0, eo), we have that the unique 
solution to ^) satisfies 

||i2(-,t) -«)(•, i)|lff.-i3xff.-i3 < Cfe"/2 

fort G [0,Toe^'^], where, 



F(x+,r) ; V b[x,t) + G{x^,t) 

Given this result, and the change of variables (H), we can immediately rewrite 
this approximation theorem in terms of the original variables. Define 

Oappix.. t) = e2 {U{e{x - t),eH) + V{e{x + t), eH)) 

+ {A{ex,et) + B{ex,et)) 

+ {F{e{x - t), e^t) + G(e(x + t), eh)) . 

Corollary 1. Fix Tq, C/ > 0, a > 13. Suppose U, V, A, B, F and G satisfy 
equations )■ Then there exists eo > and Cp > such that if the initial 

conditions for ("^ satisfy Hypothesis ^ then for e G (0,eo), the unique solution 
d{x, t) to satisfies 

\\e{x,t) ^ eapp{x,t)\\H'-^^ < CFe"/2 

/ort G [0,Toe-3]. 

Remark 2. The initial conditions for equation (jlp are obtained from those of 
simply by inverting the transformation 

The remainder of the paper is devoted to the proof of Theorem ^ and Proposition 
|l]. In the next section we give a formal derivation of equations (||)-(^. In Section 3 
we study the existence of solutions to equations (|^) and and prove Proposition 

1 . Section 4 is the technical heart of the paper and contains the proof of Theorem 

2. The proof follows the general approach for justifying modulation equations laid 
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out in ]14| , but controlling the higher order approximation requires fairly extensive 
technical modifications. In Section 5 we present the results of a variety of numerical 
computations related to Corollary 0. These computations give insight into several 
aspects of the second order approximation. First of all, it allows us to estimate 
how large the values of eo and Cp in Corollary |l| are. They also show that the 
order of e in the error estimates (i.e. 11/2) is apparently optimal. Finally, in 
the concluding section we discuss other work on second order corrections to the 
KdV approximation, both rigorous and non-rigorous, and how it relates to our own 
results. 



2. Formal Derivation of the Modulation Equations 

One can derive from (||) a system of KdV equations via the method of multiple 
time scales - this was done in [Q, for example. We extend that calculation in this 
section to include the approximating equations for the next order correction. 

To derive the modulation equations we first make the Ansatz 

u{x,t) \_of U{X^,T) \ 4^ A[X,t) + F[X_,T) 



\ v{x,t) { V{xl',T) ) [ B{x\r) + G{X~\T) ) + ^ 

where t = et, T — e^t, X — ex, X_ — X ~ t, and X^ — X + t. The two new time 
variables are the "multiple time scales" spoken of earlier. For convenience we will 
also denote u = {u,v)\ U = {U,V)\ A = (A^Bf and F = {F,Gf. 

It may seem somewhat odd that the 0{e^) correction consists of a sum of func- 
tions, as opposed to a single function. The reason for this is that for our first order 
approximation terms, U and V , we are assuming that u and v exhibit only uni- 
directional motion (right and left, respectively). A and _B, loosely, correct for the 
effect of the interaction of right and left moving waves, and they evolve on the fast 
time scale, t. There are also unidirectional second order effects, which we represent 
with F and G. Their functional form is the same as that of the first order terms 

In a moment, we will insert (^ into (||), but first we compute the effect of the 
operator A on long wavelength data. Define a function W by w{x) = W{X). We 
wish to compute \w{x), and see how that relates to W{X). The function W{X) 
represents solutions of long wavelength and hence in the Fourier domain we expect 
the frequency content of these waves to be concentrated near zero. Thus, we will 
(formally) approximate the effect of A(fc) by the first few terms of its Maclaurin 
series. 

\w{x) = |A(fc)u;(fc)} [x] 

e'^''\{k)e-~^W{k/t)dk, K = k/e 
e'^^\{eK)W{K)dK 
= :S-^[\{eK)W{K)Yx) 

= I [em + \e\iKf + le\^Kf + 0(6^)) W{K)^ (X) 

= {edx + ie^Si + ^e^ai -I- 0{e'))WiX) 
It is important to note that this approximation is only formally good to O(e^). 
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Now we insert this approximation for A and the Ansatz into (^. This is a 
necessarily messy procedure. To reduce the notation, anything formally 0(e®) or 
higher is (more or less) disregarded. Also an additional term is added to the Ansatz 
of the form: 

S\X,t) 

While this term will be treated in much the same way as the other terms in the 
Ansatz, it should be noted that this is not truly part of the next order correction. 
It will, however, be quite useful when we prove the approximation is a good one. 

We must re-express the partial derivatives in ^ in terms of the new coordinates. 
By the chain rule, we have 

dt = -edx_ + edx^ + edr + e'dr- 



Spatial derivatives of terms of the form /(X_)g(Ar_)_) or f{X)g{X±) are denoted by 
dx , though all other spatial derivatives are denoted with respect to the appropriate 
coordinate. 

So we get on the left hand side of (^) 



dt 



(10) 



u{x,t) \^3f -dx_U{X-,T) 
v(x,t) ) \ dx^V{X+,T) 

" dTU{X_,T) + drAiX,r) - dx_FiX_,T) 
dTViX+ , T) + drB{X, t) + dxM^+ ' T) 

dTF{X^,T)+drS^iX,T) 
dTG{X+,T)+drS^{X,T) 

Now we must compute the right hand side of (0). A routine calculations yields: 



RHS 



(11) 



-dx^U 
dx^V 

-^d'x.V^ 

~\dxV^ 
\dx_U^ 



-dxA 
dxB^ 



- \dx+ 



dx_F 
dx^G 

V^-dx{UV) 
dxiUV) 



--2d'xA-kd'x^F-ld%_U 



V 



-dxiUA) - dxjUF) - dxiUB) - dxiUG) 
-dx{UA) + dx_ {UF) + dx{UB) + dx{UG) 



-dx{VA) 
+dx{VA) 



So we see that we can satisfy ( 



dxiyF)^ 
dx{VF)- 

WxX 



dx{VB)-dx+{VG) 
dx{VB)+dx^{VG) 



\d\{UV) 
\dUUV) 



dxS' 
dxS^ 



(i 



formally to 0{e^ 
1 
2 



by taking 



1 



1 



OtV ^ -d-i^.V + -dx^V 
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and 

1 



1 



(i 

drB - dxB ^ ^dx_ U^{X - t, eV) + dxU{X - r, e^T)V{X + t, eV) 

Equations (^) are a pair of uncoupled Korteweg-De Vries Equations. That their 
solutions provide the first-order approximation to long wavelength solutions of ([^) 
was proven in Solutions to the KdV equations are known to exist and be 

bounded over a long time scale (see Theorem |^ above) . 

Equations (||) are a set of inhomogeneous transport equations driven by the 
solutions to the KdV equations, for which we can write down an explicit formula 
for the solutions to these equations. 

Corollary 2. The solutions to equations ^ are given by 

AiX, r)=^V'iX- r, 0)-lv'iX + r, eV) 



(12) 



whe 



(13) 



+ a{X-T, eV) + Ai(X,r) 
0) 

-/3(X + r,eV)+Bi(X,T) 



B{X, t) ^\u\X + T,G)-hj\X- T, eV) 



(14) 



Ai{X, t)^— drV^iX - T + 2s, e^s)ds 
4 Jo 

Bi{X, t) = — / drU^iX + r - 2s, e^s)ds 
4 Jo 

a{X^,T) = -e-^ I dx_ {U{X^ , s)V{X^ + 2e-'^s, s)) ds 
Jo 

T) - j ijj^x+ - 2€-^s, s)V{X+, s)) ds 

Jo 



P{X, 



Proof. The proof follows directly from Lemmas ^ and ^ which appear in the next 
section. They can also be verified by inserting the expressions for A and B back 
into (^. Furthermore, as we prove below, in spite of the prefactor of e~^, a and (3 
remain 0(1) for aU < T < Tq, for any To. □ 



The terms of O(e^) in ( ITll ) give rise to two sets of linear evolution equations, one 
for F and G and one for and S"^ . Both are inhomogeneous systems of equations 
due to the presence of terms involving t7, V, A, B and their derivatives. We 
have some freedom in the way we split up the inhomogeneous terms between these 
equations and we attempt to group them in such a way that it is easy to estimate 
the resulting solutions over the long time scales relevant for the approximation 
problem. In particular, we will break A and B up as in the explicit solutions above. 
We have: 



C. EUGENE WAYNE AND J. DOUGLAS WRIGHT 



(i ' 

where the inhomogeneous terms and are given by: 

, T) = - U{X^ , T) - (X_ , T) 



(15) 



(16) 



ax„ (^;7(x_,r)(iF2(^_^Q)_^^(^_^y)) 
i9i_«(X_,T) 



j2(X+,r) = ^9i^T/(X+,T) + i9i^l^2(X+,T) 
-i9x^l^='(X+,T) + i9i^C/2(X+,0) 



+ ax+ ( r)(^;7^(x+, 0) + T)) 



1 



(17) 



The additional terms and S*^ should satisfy 

drS^ + dxS^ —Jet + '^d + '^sp 
drS^ — dxS^ ^Jct ^ J'd ^ Jsp 



where 



J^t - -axC/(X_,eV)(G(X+,eV)+/3(X+,eV)) 

-5xC/(X_,eV) Qc/2(X+,0)- iF2(^^^^2^) 

-5x (V^(X+,6V)^^(X_,6V)) - (C/(X_,6V)y(X+,eV)) 
-9x ('y(X+, 6V)(iT/2(^_^ 0) + a(X_, eV) - i{/2(X_, eV)) 



J] = -dx{V{X+,e'T)G{X+,e\))~^d'^V\X+,e'T) + ^dxVHX+,e\) 
Jlp = -i9iv4i(X,T)-5x((;7(X_,eV) + y(X+,eV)(Ai(X,r)+i?i(X,T))) 



^dx ( eV)(ic/2(X+, 0) + eV)) 
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J2 = dxV{X+,e\){F{X^,e\)+a{X^,e\)) 
5x(X+,eV)Qt/2(X_,0)-it/2(X_,6V; 

+dx (C/(X_,eV)G(X+,eV)) + (C/(X_, eV)T/(X+, eV)) 

Jj = dx{U{X.,e\)F{X.,e'T)) + ^d'^U^X_,e'T)-^dxUHX_,e'T) 

J% = i9ii3i(X,T)+ax(C/(X_,eV) + V^(X+,eV)(Ai(X,T)+i?i(X,T))) 
+dx [u{X-,^t){-^V^{X^,{)) + eV)) 

The equations (||) are our second set of modulation equations for the terms of 
O(e^) in our long wavelength approximation. Since they are linearized, inhomo- 
geneous KdV equations, linearized about a KdV solution, they are in principle 
explicitly solvabl e |]l8| . However, the form of the solution that results is quite 
complicated (see |19| and |TT1 ) and thus it requires some effort to show that these 
solutions remain uniformly bounded in the norms which we use to bound the er- 
rors. As we noted above, the functions and 5*^ do not actually form a part of 
the approximation at 0(e*); however, we will show that they remain bounded over 
the time scales of interest as a part of controlling the error in our approximation. 



3. Estimates on the Solutions to the Modulation Equations 

Before showing that the approximation is a good one, we must first show that the 
solutions to the modulation equations are tractable in the their own right. Keeping 
in mind that our goal is show the approximation to is good for a long time, 
we need to show that solutions to the modulation equations are bounded on the 
appropriate time scale, that is, for t ^ 0{t^^). First we remark on Theorem ^ 
above. 

Notice that since T — e^t, this theorem states that we have bounded solutions 
of for t £ [0,To/e^], as we had hoped. Moreover, since the solutions to (||) 
appear in the other modulation equations (often as inhomogeneities) , that they 
are reasonably smooth and of rapid decay is crucial to showing that the other 
modulation equations are solvable over a long time, and of appropriate size. In 
particular, we will henceforth take Uq, Vq € H"'{4), where a will be suitably large. 
We now state and prove a number of lemmas. 

The first set of lemmas concerns the solutions to inhomogeneous transport equa- 
tions with zero initial conditions. From the method of charateristics, we have 
explicit formulas for solutions. 

Lemma 1. Suppose 

dru±dxu = f{X,T), m(X,0) = 0, 

with \\f{X,T)\\Hs < C for T £ [0,Toe-2]. Then \\u{-,t)\\hs < Ce'^ for t G 
[0,Toe-2]. 
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Proof. We have 



u{X,t) = 



I f{XTr±s,s) ds. 
Jo 



The integrand is bounded by C by the Sobolev embedding theorem. A naive esti- 
mate on the integral proves the result. □ 

Lemma 2. Suppose 

drU±dxu = dx.f{X±T,e^T), u(X,0) = 0. 

Then 
(18) 

u{X, r) = ±i {f{X ± r, eV) - f{X ^ r, 0)) T y ^ drfiX T r ± 2s, e's) ds 

Also, if\\f{-,T)\\H^ <C and\\dTf{:,T)\\Hs <C forTe [Q^Tq], then\\u{:,T)\\Hs < 
C/orre [0,Toe"2]. 

Proof. One can check this result explicitly. The estimate on the norm follows as in 
Lemma |l|. □ 

Lemma 3. Suppose 

dru ± dxu = 1{X + r, e^t)r{X - r, eV), u{X, 0) = 0. 
with \\1{-,T)\\hs(a) < C and ||r(-, T)||h=(4) < C for T £ [0,To], then 

u{X,t) =w(X=FT,eV) 
with |iw(-,T)||ff.(2) <C forTe [0,To] (that is for t e [0,Toe-^]). 
Proof. See appendix. □ 

Remark 3. If I and r are taken to be in H^{2), a similar proof shows that v is in 
IT^ over the long time scale. 

Remark 4. Since the proof of the Lemma does not make explicit use of the slow 
time scale dependence of the inhomogeneous factors I and r, the proof is still valid 
if the right hand is of the form 1{X + T)r{X — r, e^r), l(X + r, e^T)r{X — r) or 
1{X + T)r{X -t). 



Remark 5. A general study of the growth of solutions of the transport equation 
and related linear equations that arise in the justification of modulation equations 
was recently completely by D. Lannes p^ . 

With these results, we may now prove the estimate for A and B in Proposition 
0. That is, we have the following: 

Corollary 3. If Uq.Vq G i?°'(4), and U and V satisfy then 
sup {\\A{-,t)\\h.-^,\\B{;t)\\h.-,}<C. 

T6[0,Toe-2] 
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Proof. From Corollary ^ we know the form of A and B. By Lemma |3[ we have a 
and f3 uniformly bounded in H'^^^{2) over the long time scale. Also, by Lemma 
H we see that Ai and Bi are in the same space as OtU^ and drV^. U and V 
satisfy the KdV equations so we lose three space derivatives for the one time 
derivative here. That is, Ai and Bi are uniformly bounded in H'^^^ for the long 
time scale. □ 

We will occasionally be using an alternate, but equivalent norm, on H^{2). It is: 

s 

The associated inner product is denoted by (•, ■)h=(2) 
Lemma 4. For s > 3/2, if u E FF then 

{ud,fJ)H^<C\u\H^,\f\%.. 

Proof. The proof is similar to and simpler than that of Lemma |^, which follows. □ 
Lemma 5. For s > 3/2, if u e FF{2) then 

Proof. See appendix. □ 
Lemma 6. For f G H'{2) n 

if,dlf)H^^2)<C{\\f\\l.^,^ + \\f\\U,) 

Proof. See appendix. □ 

We may now prove the estimates on F and G in Proposition ||. That is, we have 
the following Lemma: 

Lemma 7. IfUo,Vo S i?'^(4) and U , V , A, B, F and G satisfy equations (Q/-(^, 
then F and G satisfy the estimates: 

sup {||F(.,r)||^,||G(.,T)||^}<C 

Te[0,To] 

where H = H''-^ n ^-^{2). 

Proof. The proof follows from Lemmas § - || and Gronwall's inequality. We show 
the details for F. The case for G is entirely analogous. We take the definition of 
the inner product on H to be (•, = (•, •)j?"-9(2) + Off"-^ • 

The inhomogeneity is in H"'~^{2) (the term dx_U causes the loss of deriva- 
tives) . So we take the inner product of with F and apply Lemmas |^ - 1^ and 
arrive at: 

aT||F|||<C(||F||^ + ||^^|||)<C(l + ||F|||) 
An application of Gronwall's inequality yields, for T E [0, Tq]: 

\\F\\l{T) < CTe^^ 

which concludes the proof of the Lemma and also of Proposition ^ □ 
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We now turn our eyes to the set of equations (|T^). As there are many terms 
driving these equations, many different techniques are used to show that the equa- 
tions do not blow up over the long time scale. We are aided in this task by the 
above lemmas, though certain terms will need special consideration. 

Lemma 8. Suppose U , V , A, B, F, G, and satisfy equations ^-^) and 
(p^), then and satisfy the estimates: 



sup {\\S\:T)\\H^-^o,\\S'{;r)\\H^-.o} <C 
r6[0,To£-2] 

Proof. We shall treat the equation for here. The situation for is completely 
analogous. Since equations are linear, we can consider the inhomogeneity term 
by term. First we notice that we can apply Lemma ^ to bound the growth coming 
from all terms in J^^, while Lemma ^ suffices to control all terms coming from 
J^. Thus these terms cause no growth over the long time scale. We now take a 
moment to discuss the smoothness of these terms. The least smooth term in is 
dx iV{X + r, e^T)G{X + r, e^r)). When we apply Lemma || we need to examine the 
smoothness of dT{V{X + t, £'^t)G{X + r, e^r)). Now, G is uniformly bounded in 
H'^~^, (from Lemma |^), and since G satisfies a linearized KdV equation, we have 
dTiV{X + T,t^T)G{X + T, e^r)) uniformly bounded in H"^^. However, the least 
smooth term in Jl^ is the term dx{U{X — r, e'^T)G{X + r, e^r)), which is uniformly 
bounded in iJ'^"^°(2). Thus at best is uniformly bounded in H''~^^ . 
Each term in Jj^ will require some special consideration. These terms are: 

(19) - 9x(n^ + T,eV)/3(X + T,eV)) 

(20) - dx{V{X + T,e^T)U^{X + T,Q)) 

(21) C d\A^{X,T) 

(22) G 9x(((7(^-T,eV)+T/(X + T,eV))(Ai(X,T) + Bi(X,T))) 

Terms (^9|) and (|2^) are treated with by slight variations on Lemmas || and ^. The 
technique by which ( pT] ) and (22) arc dealt with relies primarily on the prefactor of 

which appears in the definition of the functions Ai and Bi. Unfortunately, each 
computation is rather messy. 

In the case of the first of these, we apply Lemma ^ and get 

S{X, t) = -]^{V{X + T, eV)/3(X + T, eV) - V{X - r, 0)/3(X - r, 0)) 



dTV{X - r + 2s, e^s)l3{X - r + 2s, e^s)ds 



I V{X -T + 2s,e^s)dTf3{X ~T + 2s,e^s)ds 
Jq 



e 

+ Y 

The first three terms are easily bounded by the techniques discussed previously 
(namely we replace dxV with the right hand side of the KdV equation and use 
naive bounds). However, when we replace OtP, we lose the prefactor of e^. That 
is, from equation ( |32[ ) in the proof of Lemma |^, we have 

dTPiX+,T) = €-^dx{U{X+ - 2Te~^,T)V{X+,T)). 

We make this substitution into the last term of (^ to get: 

1 f 

- / V{X - r + 2s, e^s)dx{U{X - t, (i'^s)V{X - t + 2s, e^s))ds 

2 In 
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Notice that in this integral we have only terms that lie in the weighted Sobolev 
spaces, and we can use the same techniques used in the proof of Lemma ^ to 
control this term. 

The term ( pi] ) is very nearly of the form needed to apply Lemma |^. The only 
difference is that there is no dependence on e^r in one of the terms. The ideas are 
essentially the same here as in the proof of Lemma ^. Consider, 

drS + dxS = -dx{V{X + r, eV)C/2(X + t, 0)). 

The solution to this equation is given by: 

S{X,t) = - i {C/2(X + r,0)y(X + T,eV) - [/^(X - r, 0)T/(X - r, eV)} 

+ — / U'^{X-T + 2s,G)dTV{X-T + 2s,e^s)ds 
2 Jo 

If one replaces dTV{X—T+2s, e^r) in the integral by the right hand side of the KdV 
equation, and then takes naive norms, we find that this term is also controllable. 

We now turn our attention to the final two terms which involve the functions Ai 
and Bi. The calculations here are quite messy, though the ideas are straightforward. 
We replace dxV with the right hand side of the KdV equation and then apply a 
number of the same techniques used in proving Lemmas ^ and |[ The factor of 
present in the definitions of Ai and Bi is crucial. Consider 

drS + dxS 
=Cd\A^{X,T) 

=Ce^d\ [ drV^iX -T + 2s,e^s)ds 
Jo 

=Ce2a2, f (^dsidrV^iX -T + 2s,e^s))-e^dlV{X -T + 2s,e^s))ds 
Jo 

=Ce^d\ [drV^iX + t, eV) - dTV^{X - t,0)} 

+ Ce'^d\ [ d^V^iX -T + 2s,e^s)ds 
Jo 

For ease of notation, we will let P{X,t) — Ce^d\ d^V^{X — r + 2s,e^s)ds. 
Notice that by taking naive estimates on this function, we have that ||-P||hs < C 
for r € [0,Toe^^]. Thus we apply Lemma |l| to this equation to find that 5* is 
bounded on the long time interval. 

In order to deal with (^), we will rewrite dxV^ and dxU^- That is, 

drV^ = 2VdTV 

^dx(vdj,V-lidxVr + lv'^ 
= dxV. 

Where V = VdjiV-l/2{dxV)'^+2/3V^. A similar calculation yields OtW^ = dxU, 
wheieU ^ -UdjiU+l/2{dxUf-2/3U^. Notice that C/, f e i/'^-^ foj. T € [0,To], 
since they lose at most two derivatives in comparison with U and V. Similarly, we 
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have drU £ H'^ ^. So consider the equation, 

drS + dxS 

- Cdx {U{X - T, eV) + V{X + r, eV)) r) + (X, t)) 



Ce^d 



'X 



{U{X - T, eV) + V{X + r, eV)) 

- r + 2s, e^s) + U^{X + r - 2s, e^s)) d. 
dx {U{X - T, e\) + V{X + T, eV)) 
dx (V{X - T + 2s, e^s) + U{X + t - 2s, e^s)) ds 

{UiX-T, eV) + F(X + T, eV)) 



X V{X + r, eV) - V{X - t, 0) 



9t^"(^ - t + 2s, e^s)ds 



+U{X - T, eV) - C/(X + T, 0) - / drUiX + r - 2s, e^s)ds 



Notice that by taking naive estimates, the terms Q^{X,t) = dTV{X — t + 
2s, e^s) ds and Q^{X, r) = JJ" dTU{X + r — 2s, e^s) ds are uniformly bounded 
in H"^^ over the long time scale. Thus we apply Lemma |l] to the above equation 
and find that this term is well-behaved over the long time scale. □ 

4. The Validity of the Approximation 

In this section we set prove that the approximation to a true solution of made 
by the Ansatz is in fact a good one by completing the proof of Theorem ^j. 

Proof, (of Theorem ^.) 

To prove this theorem we shall need a number of lemmas. 

Lemma 9. //$ G then for e < 1, 

||A$(e.)||H. <Cei/2||$||^,+,. 

Proof. The proof here is analogous to the proof of the following lemma. □ 

Lemma 10. Let Ti{y) = y, T^iy) + l/2y^ and T^iy) ^y + l/2y^ + 3/8y^. 
Then for j = 1,3,5 if G H''+^+^ we have, for e<l, 

\\X<S>{e-)~Tj{edxMe-)\\H^ < Ce^+^Z^ 

Proof. See appendix. □ 

Now suppose that there is a solution to (||) of the form, 

(23) u{x,t) = e^^{x,t)+e^^^^R{x,t) 
where 

(24) e^^{x,t)=f^U + e^iA + F)+e'^S, 
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and R — {R^ {x , t) , Rl^ {x , t)y . We consider the term R to be the error in our 
approximation. Substituting ( p3| ) into (|^), we find that R must satisfy the equation, 

where 
(26) 

Res[e2vl/] = -9, ^^^^^ + n ^ .2,t,2 + o ^^:2,T,l ^.2 



,2\2 



y £2^2 y ^ \ )\ £2*2 y ^ 2 V A(€2*l + £2*2)2 

We have selected our modulation equations precisely so that this term is small. 
By taking the time derivative of 4" and making then making substitutions from the 
modulation equations, we find that 

Res[£2*] 

_ 2 ( {n{edx)-X)U 

,( {T,{edx)-\){A + F+\{U + Vf) 
\^{n{edx)-\){B + G+\{U + Vf) 

6 / {T^{edx) -\){{U + V){A + F + B + G) + S^) 
(27) -{Ti{edx) - \){{U + V){A + F + B + G) + S^) 

8 ( -A {2{U + V){S^ + S^) + {A + F + B + Of) 
\{2{U + V){S^ +S^) + {A + F + B + Gf) 

-\{{A + F + B + G){S^ + S^)) 
\{{A + F + B + G){S^ + S^)) 



-£ 
+£ 

+2£io 



,U -A ((^1+ 52)2) 



' A ((51 + ^2)2 

While the algebra that goes into showing this is lengthy, it should be noted that this 
step is accomplished by undoing to algebra that goes into deriving the modulation 
equations formally. 

Notice that in the above expression, all functions are of long wavelength form. 
Thus we can apply Lemmas ^ and |l^ to prove the following result. 

Lemma 11. Under the hypotheses of Theorem^, the residual satisfies the estimate: 

sup ||iZes[£2*]||^<,-i3xjj.-l3 <C£"/2 
t6[0,Toe-3] 

Notice that the loss of three more derivatives is caused by the application of 
Lemma |lO| to the term in the fourth line of equation (p7|), since , j = 1,2 are 
uniformly bounded in H'^~^^. 

We also need the following fact, 

Lemma 12. 

(i?2 - R\ A[(*1 + *2)(i^l + i?2)])^^ 
< ~ {dt{R' + i?2), (vpl + VI/2)(^1 + ^2))^^ ^ Ce'\\R\\H^^H^ 
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Proof. See appendix. 



□ 



We wish to keep the norm of R from growing too much over the long time scale. 
That is, if we can show that ||_R|| is 0(1) for t e [0, Toe"'^], we will have shown that 
our approximation is good. 

The first term on the right hand side of equation ( p5| ) will not cause any growth 
in the norm, since (/, Xf)H' = 0. The third term has the prefactor of e^^/^, which 
will assist in controlling it, and we know from the above Lemma ^ that the residual 
is small. 

If we tried to control solutions of (^5|) by applying a Gronwall type estimate to 
the time derivative of (/, f)H= , the second term would result in growth of the norm 
which would destroy our estimate over the time scale of interest. To avoid this 
problem we introduce a new energy functional which yields a norm equivalent to 
the X H'^ norm, but which does not suffer from this sort of uncontrolled growth. 

Thus we define 



That this norm is equivalent to the standard norm on H'^ x if" can be seen 
by applying the Cauchy-Schwarz inequality to the inner product, provided that we 
have e^ll^*^ + ^'■^11^= < 1. Thus we use without further comment 



We now state and prove a useful lemma. 



Lemma 13. Set s > 0. Suppose f{x),g{x) E H'' and 7(X) G H''^^ where X = ex. 
Then 



(28) EUR) - I (PI 



= X 



s + e2(i?i + R\ (^-1 + «'2)(i?i + i?2))ff.) 



7^||i?||ff^xff= < Es{R) < C\\R\\h^^hs. 



\{f{x),^{ex)g{x))H^ - {g{x),j{ex)f{x))H^\ < Ce||/|U. ||<?|U. hlk^ 



Proof. See appendix. 



□ 
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We now have all the tools needed to finish the proof of the theorem. 
= [R^ , dfR^) }jcr-i3 + [R? ^ dtR?) H"^^^ 

+ ^ {dt{R' + + Vl/2)(i?l + 

+ ^ + a^(^l + ^2)(^1 ^ ij2))^^_^^ 

< {R\dtR^)H-^^ + {R^, dtR^)H''-^^ 
+Ce^||^||//^-i3xii-<^-i3 WdtRWw-^-^xH''-^^ 

= £2 (^2 _ ^1^ ^[(,1,1 + ,I,2)(^1 ^ i?2)])^^_^^ 

+ £2 + i?2)^ (^1 ^ ^2)(^1 ^ i?2))^^_^^ 

+ £"/2 (i?2-i?l,A[(i?l+i?2)2])^^_^3 

-VCC^ \\R\\ J^(T-13 s^//tT-13 

\\dtR\\ 

+ Ce^|j^||//^-13x/i-<'-13 

+Ce^ 1 1 I ^„ _ 13 X /i-<T ~ 13 

< Ce3ij2_,3(i?) + Ce3iJ,_i3(i?) + Ce^^/^El_,,{R) 

We now state another lemma, which proves that the approximation is good over 
the long time interval. 

Lemma 14. Given C > 0, To>0, there exists eq > such that if e £ (0, Eq) and 

(29) f^{T)<C{l + 7^{T) + e''^r,^'^{T)), 77(0)^0, 
for T e [0, To] then r](T) < 2CToe^^'^'> for T £ [0, Tq] . 

Proof. See appendix. □ 

We apply Lemma |lj to the equation for the energy of the remainder and find 
that 

(30) E,^i3{R{;t))<C 

for t e [0, Toe~^]. Note that given this a priori estimate proving the existence and 
uniqueness of solutions of (psj) is a standard exercise. 
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We now use the equivalence of -Bo-- 13 to the typical norm on equation 
30|), and Lemma | to find, for t e [O^Toe-^]: 

\\u{-,t) - u;(-,t)||^<,-i3xff— 13 =\\u{-,t) ~ + e*^'5(e-,e<)||H.-i3xff— 13 

(31) =||e"/2^(-,t) + e«5(e.,et)||ff.-i3>,^.-i3 

This completes the proof. □ 

5. Some Numerics 

In this section we show the results of some numerical simulations. We performed 
these numerics to gain insight into the qualitative nature of the higher order ap- 
proximation, to estimate the values of the constants Cp and eo that appear in 
Theorem ^ and to validate the results of said theorem. 

We will choose the initial conditions of the system so that we can use the known 
solitary wave solutions to the KdV equation. We shall solve the Boussinesq equation 
numerically. Though techniques are known for finding explicit solutions to the 
linearized KdV equation (see and the resulting expressions are quite 

complicated, and so we also solve (^) numerically. 

One may wonder why we should even bother computing higher order modulation 
equations, if we have to solve them numerically. In our situation, numerically com- 
puting solutions to the Boussinesq equation is not particularly more complicated 
or time intensive than finding solutions to the linearized KdV equations. However, 
our goal is to apply these same ideas to derive corrections to the KdV approxima- 
tion for the water wave problem, whose numerical solution is a much more difficult 
task. We expect the same modulation equations to hold in these more general and 
complicated systems. Thus for the water wave problem, numerically solving the 
modulation equations should result in a great reduction in the complexity of the 
numerics. 

The solutions of (|^) and (H) are numerically computed using methods which are 
largely based around the pseudo-spectral techniques for Matlab used in . Since 
our equations are relatively simple, Matlab, though slower than other languages 
(C or Fortran, for example), performs adequately rapidly. The techniques used 
are largely built around the use of the fast Fourier transform (FFT) to compute 
the various operators and derivatives, and an iterative technique to compute the 
nonlinear terms in (|^) and the term dx{UF) in (|^). It is implicit in the time step. 

As noted previously, we use the known explicit solutions to (^. Where possible, 
we find explicit solutions for the various terms of A and B. The notable exception 
to this is in the computation of a and /3, which we compute via routine trapezoidal 
rule techniques. 

We first consider the head-on collision of two solitary waves. Note that the 
head-on collision will take place in the initial variable 9, not in either the u ot v 
variables. This is because we have (formally) decomposed the system into left and 
right moving waves when we rewrite the system as (^. We therefore take initial 
conditions such that U and V will evolve as the well-known sech-squared solitary 
wave solutions to (||). That is we take: 

u{x,0) =6e2sech2(ea;-10) 

v{x,0) =6e2sech^(ea; + 10) 
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Figure 1. Initial profile for head-on collision, e = 0.1 
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Figure 2. Head-on collision, e = 0.1 

as initial conditions. 

Figures |l| and || show the solution to the Boussinesq equation, as well as the 
KdV approximation and the second order correction on the same plot, at the start 
and at the collision. Here e = 0.1. We remark on several features of the Boussinesq 
equation that are not reflected in the KdV approximation, but are present in the 
second order correction. 

First, in the KdV approximation, during the collision, the two waves add in linear 
superposition (this can be seen as the the KdV equations evolve independently). 
However, the solutions to the Boussinesq equation do not display this simple linear 
property during the collision; the total height of the wave is slightly less than 
the sum of the two heights of the two waves independently. The second order 
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Figure 3. Close-up of peak of waves, during the head-on coUision. 
e = 0.1 
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Figure 4. A "shadow wave" and dispersive wave train in the 
head-on coUision. e = 0.1 

correction does a notably better job at displaying this feature (see Figure ||). The 
second feature we notice is the presence of "shadow waves" with dispersive wave 
trains (see Figure ^) in the solution to (|l]). These are not present in the KdV 
approximation, but are seen in the second order correction. 

From these pictures, we see that the second order correction is in fact doing a 
better job than simply the KdV approximation alone. In order to quantify this, we 
computed the solution for variety of values of e, and computed the value of the 
and L°° error of the KdV and second order approximations. The time to collision 
is of 0(e~^), and on this time scale and slightly beyond, the maximum error occurs 
during the collision. 
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Figure 6. sup \\u — w\\l=° vs. e for head-on collision. 









KdV 


3.49 


3.93 


KdV + second order correction 


5.55 


6.04 



Table 1. Order of the approximation, numerically computed, for 
the head-on collision 



Figures |^ and || display log- log plots of the and L°° error versus e respectively. 
The slopes of these lines are the order of the correction. We note that we have 
used only those values of e < 0.1 in computing these slopes, as we expect the error 
estimates to hold if e is sufficiently small. Table B summarizes the results. 
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KdV 


35.5 


25.8 


KdV + second order correction 


503 


523 



Table 2. Value of Cf, numerically computed, for the head-on collision 




Figure 7. Overtaking Wave before the CoUision, e — 0.05 



From this we see that our estimate of the error made in approximating the true 
solution by the second order approximation is optimal, in terms of powers of e. By 
taking note of the y-intercept of these lines, we can get an estimate on the value 
of the constant Cf in each case, see Table (j^). Unfortunately, these values of the 
constant are quite large for the second order correction. We also note that it is not 
so much the actual value of the leading coefficient that matters, as is the location 
(in e) at which the second order correction and the KdV correction return the same 
error. That is, graphically, where the lines in figures (||) and (^) cross. 

The next simulation was that of right moving overtaking waves. We take initial 
data such that U will evolve as the famous two soliton solution to (|4|). Since we 
are not interested in left moving waves, we take initial data for v to be zero. Note 
that V does not remain zero, however, due to the coupling. 

Unlike the previous situation, the time scale of the overtaking wave collision is 
0{e~^). To observe the entirety of the collision, we take Tq = 8. We also observe 
that the error in the approximation is largest at the end of the interval [0,Toe^^]. 
From the proof of Lemma one can see that as Tq increases, eo decreases. This 
requires smaller values of e, which in turn necessitates running the simulation for a 
longer period of time. 

Figures |^ and ^ display the values of u and the approximations at various 
times during the collision. As in the case of the head-on collision, the second order 
correction picks up the presence of a dispersive wave, which is not seen in the KdV 
approximation (see Figure |o|) 

It is well known that in the two soliton interaction, the waves are phase-shifted 
after the collision (that is, the faster wave is further ahead after the collision than it 
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Figure 8. Overtaking Wave at the Collision, e = 0.05 




300 350 400 450 

X 



Figure 9. Overtaking Wave after the Colhsion, e — 0.05 



would have been had no interaction taken place, and the slower wave falls behind in 
a similar fashion). Overtaking waves in the Boussinesq equation share this feature, 
though with a different phase shift. This can be seen in Figure ^, where the KdV 
approximation is leading the Boussinesq solution. The second order correction 
noticeably "fixes" this problem. In Figures |ll] we plot the locations of the peaks. 
Note that these figures reflect the fact that the numerics are computed in a moving 
reference frame (moving to the right with unit velocity). 

In Figure |l^ we plot the error in the phase shifts for the two approximations 
versus e. Notice that the slope for the second order correction is steeper than that 
of the KdV approximation 
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Figure 10. Dispersive Wave for the Overtaking Wave 
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Figure 11. Wave Peak Locations, e — 0.05 









KdV 


3.43 


3.88 


KdV + second order correction 


5.42 


5.86 



Table 3. Order of the approximation, numericaUy computed, for 
the overtaking colhsion 



In Figures ^ and ^ we plot the maximum of the and L°° error for the two 
approximations versus e on a log-log plot (as we did for the head-on interaction 
earlier). We summarize the results in Tables (||) and (^). 
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Figure 12. Error in Phase Shift vs. e 
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Figure 13. sup \\u — w\\i^2 vs. e for Overtaking Wave Collision 







L°° 


KdV 


1300 


860 


KdV + second order correction 


19,500 


19, 900 



Table 4. Value of Cf, numerically computed, for the overtaking collision 

6. Conclusions 

We conclude by briefly surveying other work on the derivation of higher order 

modulation equations for water waves and related systems. 

For the actual water wave equations there have been a number of studies of 
corrections to the KdV approximation to water waves spanning the spectrum from 
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Figure 14. max ||u — w||loo vs. e for Overtaking Wave Collision 

non-rigorous asymptotic expansions ||] , [Q , to numerical solutions of the equa- 
tions of motion and comparison with the KdV predictions , , |^ to experimen- 
tal investigations jl6| , [|| . We concentrate here on the theoretical studies since they 
have the closest connection to our work. In the investigations of Byatt-Smith 
11 and Su & Marie || the focus is on the head-on collision of solitary waves. This 
has several consequences. First of all the authors assume that the initial conditions 
are of a special form, namely a pair of counter-propagating solitary waves. The 
higher-order corrections to the solution then exploit this special form by including 
not only a correction to the amplitude of the solution, but a phase shift for each 
wave as it undergoes the collision. This is a very reasonable hypothesis in these 
physical circumstances, but one which can't easily be adapted to the more general 
type of initial conditions considered in our work. Furthermore, since these papers 
consider specifically the head-on collision of solitary waves they are concerned with 
events which occur on relatively short time scales (i.e. time scales of O(^) in our 
scaling.) As noted in (Q, p. 503) these expansions are not uniformly valid in 
time and it is not clear whether or not their solutions could be controlled over time 
scales of 0{-^). It is worth noting that in spite of the differences between our ap- 
proach and those discussed here, Byatt-Smith Q| also finds that corrections to the 
amplitude of the solitary wave evolve according to the linearized KdV equation. 

Another set of papers by Sachs, Q, Zho and Su [0, and Haragu§-Courcelle, 
Sattinger, and NichoUs considers corrections to the KdV approximation for 
unidirectional motion. The first two of these papers study this question in the 
context of water waves, while [ [lO| | studies the KdV approximation to solutions of 
the Euler-Poisson equations. The focus of these papers (particularly [|9) and [|To|) 
is rather different than ours, however. Both derive an inhomogeneous linearized 
KdV equation for the correction to the KdV approximation. However, rather than 
deriving rigorous estimates of the difference between the approximate solutions pro- 
vided by the model equations and the true solutions they focus on the nature of the 
solutions of the linearized, inhomogeneous KdV equation. In particular, Sachs [ p^ 
shows that if one linearizes about the A^-soliton solution of the KdV equation the 
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resulting inhomogeneous equation has solutions which have no secular growth. In 
|To| the authors obtain explicit solutions of the linearized KdV equation, particu- 
larly for the case in which one linearizes about the two-soliton solution of the KdV 
equation. Haragu§-Courcelle, et al. then compare the approximation they obtain 
to numerically computed solutions of the Euler-Poisson equation and they find that 
the addition of the solution of the linearized KdV equation to the approximation 
given by the two-soliton solution of the KdV equation does a significantly better 
job of approximating the solution of the Euler-Poisson equation. In particular, they 
note that the prediction of the phase shift that occurs when a "fast" traveling wave 
overtakes a slower one is significantly better when the second order correction is 
included. This effect is also present in our approximation - see Figure (p^. 
Acknowledgments: The work of the authors was supported in part by the NSF 
grant DMS-0103915. The authors thank D. NichoUs and D. Sattinger for several 
discussions of their work prior to publication. The second author also thanks G. 
Schneider for numerous discussions of the derivation and justification of modulation 
equations and H. Segur for discussions of higher order approximations to water 
waves. 



Proof. For Lemma |3|: We consider the case with the "minus" sign on the left 
hand sign for simplicity. The other case is analogous. First we change variables 
to = X + T, T = e^T, and v{X+,T) = u{X, r). Under this change we get the 
equation: 



7. Appendix 



(32) dTv{X+,T) = e-^l{X+,T)r{X+ ~ 2Te~\T) 

This can solved by integrating with respect to the variable T. We get. 



v{X+,T) = e-^ f 1{X+, s)r{X+ - 2se-^2, s)ds 




<e-^ [ {l + Xl)\l{X+,s)\{l + {X+-2se-YMX+-2se-\s)\ ds 




ds 



ds 




< C||?||H-,(4)||r||^.(4)arctan(2re-2) 
<C\\l\\HHi)\\r\\HHi) 



□ 
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Proof. For Lemma ^ 

{ud,f,)HH2) = ("5./, /)ff.-i(2) + ((1 + X')d',{ud,f), (1 + X^)d'J)L2 



<c\u\Hs^2)\f\'Hs^,^ + - j{i + x'YudAdUYdx 

<C\u\Hsi2)\f?Hs^^)-\ J dMl + xYuWjfdx 

<C\u\Hs^2)\f\h^,^-l J dM{l + ^^)diffdx 
xu{l + x^){dlffdx 

<C|U|^.(2)|/|1^.(2) 



Proof. For Lemma ^ In this proof we use the standard norm on H^{2). 
{f,dlf)H^(2) - -6(a:2/,a,/)^. - 6((1 + 

< C^(ll/ll//=(2)ll/lk^" + + C{\\f\U^(2)\\xdlf\\H^) 

So now consider 

|2fl|2 _|i c,2f||2 , II c)S-1/„a2 J.MI2 , 11 £)S c,2 j-mi 2 



□ 



\xdtfrH. = 11x9^/11^.-. + \\dl-\xdtf)\\i. + WAxdtf)\ 



L2 



We now treat the last term in the above, as the middle term can be handled in a 
similar fashion, and the first is easily dealt with. 

\mxdlf)\\l. <C\\dl+'f\\l. + 

<C||/|ll^.+4+ j x'd^+'fd^^+'fdx 

<C||/||2,.+, + J x^dl+^fdlfdx + A j xdl+^fdlfdx 

+ 2 j dl+^fdlfdx 

<C{\\f\\ls+. + \\f\\H^^A\f\\HH2)) 

This estimate completes the proof. □ 



Proof. For Lemma 10: Notice that the polynomials Tj are the first, third and fifth 
order polynomial expansions of y/ \J^~ about y = 0. Moreover, note that only 
odd powers appear in the expansion. So, by Taylor's theorem, there is a constant 
C such that |y/yr^- T,(y)| < 

We shall now use the Fourier transform version of the Sobolev norms in the 
following computation, which concludes the proof. Consider 
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\me-)-T,{edxMe-)\\%. 

(1 + k^lX^S) - Tj{edxj^{ex)\^dk 



ik 



-Tj{ik))^k/e)\^dk 



< Ce-^ J{l + ky\k^+^^{k/e)\^dk, K = k/e 
= Ce^^+^ J (1 + {eKfy\K^+^^{K)\^dK 

< Ce^^+^ j {\ + Ky\K^+^^{K)\'^dK 



□ 



Proof. For Lemma 12 



-11/2 



(Res[*]i + Res[*]2, + ^'2)(^i 4. R^)) 



□ 
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Proof. For Lemma 13: We shall be using the fact that, by the Sobolev embedding 
theorem, we have 7 and its first s derivatives in L°°. 

\{f{x),l{ex)g{x))H' ~ igix),j{ex)f{x))H^\ 

S 

=\Y.{dif{x),di{^{ex)g{x)))^, - {dig{x),di{j{ex)fix)))^, \ 
j=o 

^lEE^J' {idifi^),dUi{exm-'gix)))^, - {digix),dUliexm-' f{x))) | 
j=o 1=0 

= lEEc,7 {dif{x)di-^gix)-d^gix)di-'fix),dUl{ex)))^, \ 
j=o 1=0 

=e|EEc,7 {difix)di-'gix)-d=gix)di-'fix),e'-'d'^jiex))^,\ 
j=o 1=1 

s j 



EE^^' / {dif{x)di-^g{x)-d^g{x)di-^f{x)) {e^-^d'^j{ex)) dx\ 
j=o 1=1 

<e||7lk».~ EE^^-' / I idif{x)di-'g{x) - g{x)di-' f{x)) \dx 
<CeMw^M\f\\HA\9W 

□ 



Proof. For Lemma 14: Functions which obey the inequality are bounded above 
by solutions to the following family of ordinary differential equations, 

r)(r; e) = C{1 + r^T; e) + e-'^rf^'^iT- e)), ,7(0; e) = 0, 

and so we prove the result for these equations. 

By separation of variables, we have that rj{T] 0) — e'-^^ — 1. We notice that for 
fixed T, TjiT] e) is a continuous and increasing function of e. This follows since 
solutions of ODEs depend smoothly on their parameters and that the right hand 
side of the differential equations is increasing in e. 

Thus, by the intermediate value theorem, there exists eg such that ri{TQ]tQ) = 
e^^ . Moreover since is a decreasing function for e > 0, we have rjiTo] e) < e~'^ 
for e e (0,eo). We further note that for fixed e, 77 is continuous and increasing in 
T. So we have r]{T- e) < e^^ for T G [0, To] and e e (0, eo). 

Thus we have, 

77(r;e)<C(l + 2r;(r;e)) r;(0;e) = 0, 

for T E [Q,Tq] and e E (0,eo). We apply Gronwall's inequality to this to prove the 
result. □ 
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